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Abstract 

We study the regularization dependence on the quenched Schwinger-Dyson equations in general 
gauge by applying the two types of regularizations, the four and three dimensional momentum 
cutoffs. The obtained results indicate that the solutions are not drastically affected by the choice 
of two different cutoff prescriptions. We then think that both the regularizations can nicely be 
adopted in the analyses for the Schwinger-Dyson equations. 
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I. INTRODUCTION 


The chiral symmetry breaking in the strongly coupling system is interesting phenomena 
in quantum chromodynamics. Although the coupling strength is small in realistic quantum 
electrodynamics (QED), the chiral symmetry can be broken when we consider strongly 
coupled QED. Then it is interesting to study such a system in Abelian gauge theories. 


For the investigation on the above mentioned chiral symmetry breaking coming from 


Te analysis based on the Schwinger-Dyson equation (SDE) may 
l|. The SDE is the set of equations for green functions whose 
brmation on the field renormalization and dynamical mass generation 
]). The analysis of the SDE should conceptually lead the gauge and 
regularization independent solutions, since the equations stem from a renormalizable gauge 
theory. However, the solutions practically depend on the chosen gauge, because the equations 


dynamical mass generation, t 
be an appropriate approach 
solutions give the in 
(for reviews, see, [2 


are derived through applying some approximations in the intermediate steps. The extensive 
analyses on the gauge dependence are given in Q where the effects of generalised vertices are 
studied in unquenched QED with four dimensional cutoff method. Concerning on the similar 


genera 

e,„ B 


isations of the equations, a lot of works have been done with various approaches, see, 


7], A 


procedures [§], since the quenched SDE can be regarded as the generalised version of the 
Nambu Jona Lasinio-type gap equations in which regularization has effects on the model 
predictions [ 9 ]. We then think it may be interesting to study the regularization dependence 
on the solutions of the SDE. 


so, it is known that the physical predictions as well depend on regularization 


In this letter, we shall numerically solve the SDE in general gauge with two regularization 
procedures, the four dimensional (4D) and three dimensional (3D) cutoff regularizations, 
then make the comparison between these two methods. The importance of studying the 
equations with the 3D cutoff regularization lies on the fact that the solutions can smoothly 
be continued to the ones obtained in the finite temperature system, since the equations are 
usually investigated by using the 3D cutoff scheme at finite temperature 10]. 


This paper is organised as follows; Section |TT] presents two types of equations. We show 
the numerical results for the field renormalization factor and the dynamically generated 
mass in Sec. IIIII The concluding remarks are given in Sec. IIV1 
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II. SCHWINGER-DYSON EQUATION 

The Schwinger-Dyson equation for fermion self-energy E(P) is written by 

S(P) = ie 2 J ^7 “D^(P - Q)S(Q)r (P, Q), ( 1 ) 

where e is the coupling strength, D^ U (P — Q ) and S'(Q) are the gauge boson and fermion 
propagators, and T U (P,Q) is the vertex function on the gauge boson and the fermion. In 
this letter, we use the capital letters as P, Q for expressing the four dimensional momenta, 
namely, P M = (p 0 ,p) and Q lt = (go,q)- For D^ w and T", we employ the following tree-level 
forms 

D„„(K) = ^f + ( ( 2 ) 

r-(p,Q) = Y, ( 3 ) 

with the gauge parameter £ and K tl = P M — Q^. 


A. Equations with the four dimensional cutoff 


In the SDE with four dimensional cutoff, we define the fermion propagator by 

S(Q) = 1 


_ (4) 

A(Q)Q,r-B(QY 

where A(Q) and B{Q ) indicate the field strength factor and the mass function. The insertion 
of these quantities leads 


A(P)P^ - B(P ) = 
P* 7 m - + e 2 


d A Q 


IT 


+ (1 - £) 


A{Q)Q P r ~ B(Q) 


7 


(5) 


K 2 v 10 

with the fermion bare mass mo appearing in the Lagrangian. Taking the trace after multi 
plying P p 7 p and without the multiplication give the equations for A(P) and P(P) as 


MP) = 1- T2 

P(P) = m 0 — e 


d 4 Q 

i( 27r) 4 

2 I d 4 Q 


, P'Q , ( [P-K)(Q K ) 

S+ — 


7 *(2vr) 4 
where = 1 ± f, £3 = 3 + £ and 


6 


10 
1 ' 
10 


10 


A'(Q)A(Q), 


A '(Q)B(Q), 


A '(Q) = 


-1 


A 2 (Q)Q 2 -B 2 (Q)' 


( 6 ) 

(7) 

( 8 ) 
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Performing the angular integration after the Wick rotation, we obtain the following equa¬ 
tions, 


A(P) 
B{P ) 


1 + 


m 0 


47T 


(■Air 


'<5 2 

°4D 


d Q 2 


47T 


r-A 2 


P 4 

r^2 


dQ 2 


'5!d 


0(P - Q) + 0(Q - P) 


A(Q)A(Q), 


P 2 


0(P - Q) + 0(Q - P) 


A(Q)fl(Q), 


(9) 

( 10 ) 


with 


A (Q) 


1 

A 2 (Q)Q 2 + B 2 (Q)' 


( 11 ) 


where a = e 2 /(47r) and we introduce the ultraviolet and infrared cutoffs, A 4D and <5 4 d- These 
are the equations with the four dimensional cutoff scheme. 


B. Equations with the three dimensional cutoff 


In the three dimensional cutoff, we need to consider the following fermion propagator 


S(po,p) = 


C(Po,p)loP° + A(p 0 ,p)j i p i - B(p 0 ,p) 


( 12 ) 


with p — |p|, since the time and space directions should be treated separately. After the 
Wick rotation and a bit of algebras one obtains the following forms 


C(p 0 ,p) = 1 + 

A(po,p) = 1 + 
B(p 0 ,p) = rri 0 


a 

27T 2 , 
a 

2n 2 , 
a 


2n 2 


j-OO rA 3D 

/ dq 0 dq [l C AA(q 0 ,q)+lccC{q 0 ,q)]A(q 0 ,q), 

OO J &3T) 

roo rA 3 D 

I dqo dq [l AA A(q 0 ,q)+l AC C(q 0 ,q)}A(q 0 ,q), 

—oo 4<5 3D 

poo !■ A 3D 

dqo dq [l B B(q 0 ,q)]A(q 0 ,q), 
o J <5 3 e> 


(13) 

(14) 

(15) 


with 


A(q 0 ,q) = 


C 2 (qo, q)q'o + A 2 (q 0 , q)q 2 + B 2 (q 0 , q ) 


(16) 
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and 



( 17 ) 


(19) 


(18) 


( 20 ) 


Zb = &7i, 

/, _ 1_ i n + (9 ~ p) 2 


( 21 ) 


2 p kl + (5 + p ) 2 ’ 


( 22 ) 



2p L*S+ (9 - p) 2 ^0 + (9 + p) 2 


(23) 


where fco = po — Q'o- Thus we need to consider the three unknown functions C, A and B with 


two variables p 0 and p, then the numerical analyses become much more difficult comparing 
to the above mentioned four dimensional cutoff case. 

III. NUMERICAL SOLUTION 

In this section, we numerically solve the equations with the four and three dimensional 
cutoff procedures by using the iteration method, then make the comparison on the obtained 
results. 

A. Solutions with the 4D cutoff scheme 

We show the numerical results of A(P 2 ) and B(P 2 ) in Fig. [0 One sees that A(p 2 ) 
increases with respect to £. This can easily be understood because A has the form of 
A = 1 + 0(£), so it becomes larger when £ increases. On the other hand, B{p 2 ) decreases 
when £ becomes larger. This can also be understood by following discussion; B has the 
form of B oc (3 + £) J dQF(Q)A(Q), and although 3 + £ becomes larger with increasing £, 
A (Q) = 1/(A 2 Q 2 + B 2 ) decreases when A increases. Consequently, B has smaller value for 
larger £. Note that A = 1 always persists in the case £ = 0 as obviously read from Eq. flS]), 
which is the well-known consequence of choosing the Landau gauge. We also studied different 
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FIG. 1. Gauge dependence for a = 2.5 and 3.5, with mo = 0, $ 4 d = O.OIA 4 D. 

values of a, and found that the above mentioned tendency did not change, so we only showed 
the results with a = 2.5 and 3.5 here. 

It may also be worth studying the case with finite m 0 . Figure [2] shows the numerical 
results with mo = O.IA 4 D. One notes that the obtained values of A and B are closer for 




FIG. 2. Results for a = 2.5, with m 0 = O.IA4D, &id = O.OIA4D. 

various gauge comparing to the case with m 0 = 0. We can read from the results that the 
value of B is dominated by the factor 3 + £ rather than A (Q) because A has closer values 
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for this case as seen in the left panel of Fig[2l 


B. Solutions with the 3D cutoff scheme 

Here we show the numerical results of C(po,p), A(p 0 ,p ) and B(po,p) for various values 
of £ with the three dimensional cutoff. Figures [3] and [4] display the solutions for a = 2.5 




FIG. 3. Gauge dependence for a = 2.5, with mo = 0 and do = 0.01A 3 d. 

and 3.5 with the gauges, £ = 0, 1 and 3. In performing the integral in go direction, we set 
the lower and upper limit as f A ^ o dqo with A 0 = 10A 3 d because it is technically difficult to 
directly take the infinite range. We have numerically confirmed that the solutions are not 
affected by the choice of the cutoff A 0 if we take large enough value for it such as A 0 > 5A 3D . 

From the obtained curves, one sees the similar tendencies that the renormalization factors 
A and C are large when £ is large, while the mass factor B is small for larger £. It may be 
interesting to note that A and C are close to 1, but deviate from 1 in the case of the Landau 
gauge, £ = 0, which comes from the effect of the separation of the 4D momentum P to 3D 
momentum ( po , p). 

We next check the effect of the finite bare mass m 0 7 ^ 0. Figure [5] shows the solutions 
with mo = 0.1 A 3 d and a = 2.5 for various gauges. Again, one sees the similar qualitative 
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FIG. 4. Gauge dependence for a = 3.5, with mo = 0 and 63 = O.OIA 3 D. 



FIG. 5. Results for a = 2.5, with mo = O.IA 3 D and £3 = O.OIA 3 D. 


tendency with the 4D case; the renormalization factors A and C become close to 1 comparing 























to the results with m® = 0, and the mass factor B increases with £. 


C. Comparison between 4D and 3D cutoff regularizations 


We have seen above that the qualitative feature of the solutions with the 4D and 3D 
cutoff regularizations are similar. It may also be interesting to show the quantitative com¬ 
parison between these regularizations, although the direct comparison is not possible since 
the variables are different in two methods. In Figs. [6] and [71 we align the results of A(P) 
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FIG. 6. A(P), B{P) in the 4D cut (left), and A(0,p), B(0,p) in the 3D cut (right) for a = 2.5. 


and B(P) as the functions of P in the 4D cutoff, and A(0,p) and B(0,p) as the functions 
of p in the 3D cutoff. 

We find that the results between two regularizations do not alter considerably, while the 
deviations between different gauges are rather large. Then we read that, concerning on the 
solutions on the momentum direction, the regularization dependence is not serious compare 
to the gauge dependence. 
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FIG. 7. A(P), B{P) in the 4D cut (left), and A(0,p), B(0,p) in the 3D cut (right) for a = 3.5. 


IV. CONCLUDING REMARKS 


We have studied the regularization dependence on the quenched SDE in general gauge 
through applying the four and three dimensional cutoff methods in this paper. The charac¬ 
teristic technical difference lies on the number of variables where only one momentum, P, 
exists in the 4D cutoff, on the other hand there appears two different variables, p 0 and p, 
which makes the numerical analysis challenging. The obtained results between the 4D and 
3D cutoff methods show that the regularization dependence on the solutions is not drastic 
comparing to the one on the gauge parameter. This indicates that both the regularization 
prescriptions can nicely be adopted for the analysis on the SDE. 

The gauge dependence seen in the results are due to the applied approximations in deriv¬ 
ing the equations. Especially, the tree-level approximated form of the photon propagator in 
Eq. (|2]) is crucial when we consider the gauge dependence since the gauge parameter mani¬ 
festly appears in the equations. Therefore, for the sake of obtaining the gauge independent 
solutions as indicated by gauge theories, the general analyses, such as the ones based on the 
unquenched equations [5, 6], and the generalized vertices jzj], are important. 
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